Optical addressing of an individual erbium ion in silicon 
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The detection of electron spins associated with single defects in 
solids is a critical operation for a range of quantum information and 
measurement applications currently under development ^. To date, 
it has only been accomplished for two centres in crystalline solids: 
phosphorus in silicon using electrical readout based on a single elec- 
tron transistor (SET)^ and nitrogen-vacancy centres in diamond us- 
ing optical readouP^. A spin readout fidelity of about 90% has been 
demonstrated with both electrical readout^ and optical readouP^E3 
however, the thermal limitations of the electrical readout and the poor 
photon collection efficiency of the optical readout hinder achieving 
the high fidelity required for quantum information applications. Here 
we demonstrate a hybrid approach using optical excitation to change 
the charge state of the defect centre in a silicon-based SET, condi- 
tional on its spin state, and then detecting this change electrically. 
The optical frequency addressing in high spectral resolution conquers 
the thermal broadening limitation of the previous electrical readout 
and charge sensing avoids the difficulties of efficient photon collec- 
tion. This is done with erbium in silicon and has the potential to en- 
able new architectures for quantum information processing devices 
and to dramatically increase the range of defect centres that can be 
exploited. Further, the efficient electrical detection of the optical ex- 
citation of single sites in silicon is a major step in developing an 
interconnect between silicon and optical based quantum computing 
technologies. 

The potential for a hybrid optical/electrical single- spin readout 
was recently established by Steger et al., revealing long nuclear spin 
coherence time with an ensemble of P ions in highly purified ^^SM 
The readout of the spin ensemble was demonstrated by detecting the 
photocurrent generated when the excitonic transition at 1,078 nm as- 
sociated with the P ions was excited. In the current work we demon- 
strate single- site detection by electrically detecting the optical exci- 
tation of the ^Ii5/2 - ^Ii3/2 transition of single Er ions implanted into 
silicon, resolving both electronic Zeeman and hyperfine structure. 
The efficient readout required for single- site detection is achieved by 
measuring the photo-induced change in the site's charge state using 
a SET, rather than detecting the associated photocurrent. 

Erbium's large electronic magnetic moment of its ground state, 
the 7=7/2 nuclear spin of its 167 isotope along with the coincidence 



detection efficiency and low frequency resolution making state read- 
out not feasible. 
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of silica optical fibers make Er centres appealing for quantum infor- 
mation applications^^ ^^. The few existing studies in samples with 
a high concentration of Er have shown 0.1s nuclear spin relaxation 
time^^ and 100 jis electron spin dephasing time^^. Single-site detec- 
tion grants access to low Er densities where one expects drastically 
enhanced coherence times in analogy to the impressive recent de- 
velopment for P in Si^^. The low emission rate from optically ex- 
cited rare-earth ions s uch a s Er^^ makes pure optical detection of 
single sites challenging^ffi! Recently though, the optical detection 
of a single rare-earth ion was demonstrated in a Pr-doped YAG nano- 
crystal^^. The technique employed, involved the 2-step excitation of 
the ion to a high lying 5d-electron state and detecting the resultant 
emission^^. It was conducted at room temperature and exhibited low 



Ii3/2 transition is between states within the inner 4f- 
electron shell of the Er^+ ion, which is well shielded from the sur- 
rounding lattice by filled outer shells, resulting in narrow spectral 
linewidths and the potential for high resolution frequency addressing. 
At liquid helium temperatures homogeneous linewidths as narrow as 
50 Hz have been observed for the transition in Er^^:Y2Si05^^. Prior 
to the present work there have not been any sub-inhomogeneous 
studies conducted on optical transition in Er centres in silicon. The 
observed lifetime of emission of 2 ms from the "^113/2 state for Er^^ 
ions in silicon implies a minimum linewidth of 150 Hz^^. 

The resonant photoionization of individual Er^^ ions is studied in 
an Er-implanted SET (Fig. [T^), which works as a charge sensor. The 
"^115/2 - "^113/2 transition of an Er^"^ ion has a relatively high probabil- 
ity, when a laser is tuned to its resonant wavelength, and the Er^^ ion 
could be further ionized due to a second-photon process or an Auger 
process. The charge displacement induced by an ionization event si- 
multaneously leads to a change in the tunnelling current of the SET. 
To get a high sensitivity, the SET is biased close to the degeneracy 
point between two charge states, i.e. at the edge of one Coulomb 
peak (Fig. [TJ)). Accordingly, the transconductance is large, and a 
small charge displacement in the sensitive region will lead to a sig- 
nificant change in the tunnelling currenl^^^^ The photoionization of 
individual Er^^ ions leads to a significant change in tunnelling cur- 
rent (Fig. flp). The "^115/2 - "^1^/2 transition of each Er^+ ion has a 
specific resonant photon energy, so individual Er^+ ions are distin- 
guished by the resonant photon energy. When the laser is tuned to a 
non-resonant wavelength, the tunnelling current mainly stays at the 
background level as shown in Fig. IT]:. In contrast, when the laser is 
tuned to a resonant wavelength of an Er^+ ion, the photoionization 
of the Er^+ ion leads to a rise in the tunnelling current, and then the 
current drops back due to its neutralization, contributing a two-level 
current- time trace (Fig. [iji), which suggests only one single Er^^ 
ion is ionized (details in Methods). Figure [it shows a photoioniza- 
tion spectrum of a single Er^+ ion. Current-time traces are recorded 
at a series of photon energies, and then the histogram showing the 
distribution of current in time is plotted as a function of the photon 
energy detuning. The colour in Fig. [T^ represents the time (J^tbin) 
during which the current stays within one bin, and a 0.02-nA bin size 
is used for all the analysis. 

As shown in Fig. [T^, the SET has a Si channel wrapped with 
the gate. The SET is biased below the threshold voltage, so that the 
current tunnels through the corner regions of the Si channel^^. Con- 
sequently, the charge sensor is more sensitive to the Er^^ ions which 
are closer to the corner regions in the channel, and different Er^^ ions 
have diff'erent capacitive coupling leading to diff'erent detection sen- 
sitivity. The change in current (Fig. [TJD-d) accords with the loss of 
an electron, indicating that it's due to the ionization of the Er centre, 
whereas the gain of an electron will lead to a shift opposite to that 
in Fig. fTt. The small fluctuations in current, which we attribute to 
the trap states in the insulating layer or the oxide layer with weak ca- 
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FIG. 1. Photoionization spectroscopy of an individual Er^+ ion. a, Coloured scanning electron micrograph of a typical SET device used in 
this study and a band structure of Er^+ ions in silicon. Top blow-up: Schematic cross-section of the SET showing the optical addressing of 
individual Er^^ ions, b, The SET charge sensing scheme. The loss of an electron due to photoionization induces a transient shift of the I-Vg 
curve towards lower gate voltage, causing a change in current from I(q^) to I(q^). c, d, The current-time traces recorded with a fixed gate 
voltage (Vm) under (c) non-resonant and (d) resonant illumination, e, The histogram of current-time traces as a function of the photon energy 
detuning. The photon energy of the illumination is detuned with respect to the centred wavelength of 1,537.9 nm. 



pacitive coupling^^, can be suppressed by a proper anneal before the 
device fabrication. The readout efficiency is mainly limited by the 
efficiency of the excitation from the "^113/2 excited state into the con- 
duction band, which can be increased to close to 100%, by increasing 
the intensity of the light used to drive this final ionization step. We 
observe resonances via the photoionization spectroscopy mostly be- 
tween 1,535 nm and 1,539 nm, which is consistent with the ^Ii5/2 - 
^Ii3/2 transition of Er^^ ions in silicon^^ ^^. 

In the next experiment, we study the Zeeman effect of individ- 
ual Er^+ ions, as the Zeeman eff'ect is an essential tool to determine 
the site symmetry of Er centres. Er^^ ions tend to take 3+ valence 
characteristic of the Si lattice, so the 4f electrons of Er^+ ions have 



and the first excited state of ^Ii3/2^^- 
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the ground state of "^115/2 
degeneracy is lifted by the crystal field, so that each state splits into 
several levels depending on the symmetry of the Er centre^^. The 
transition between the lowest level of ^Ii5/2 and the lowest level of 
"^113/2 is responsible for the strong emission band around 1.54 jim, 
and the Zeeman splitting of those two levels in the case of double 
degeneracy is shown in Fig. [2^. The doublet states can be described 
by an eff'ective spin 5 = 1/2, and the Zeeman interaction has the form: 
H-Pe^ • g • S, where pe is the electronic Bohr magneton, B is the 
magnetic field, and g is the g-f actor matrix^^. The Zeeman splitting 
energy of the higher (lower) energy doublet is proportional to gu 
(gi). As shown in Fig. [2^, the energy diff'erence between the two 
AMs=±l (AMs-0) transitions can be described by AE - PeAgB, 
where Ag is the g-factor difference. In this study, we measured the 
Zeeman splitting of 4 spectrally isolated Er resonances, and observed 
the g-factor difference from 1.6 to 10.8. 

Figure 03,c shows the Zeeman splitting of Er^^ ions. Current- 
time traces are taken at a series of photon energies and magnetic 
fields, and each pixel in Fig. [2]3,c represents one current-time trace. 
When an Er^+ ion is ionized, the current will exceed a certain thresh- 
old, which is determined by the background current fluctuation un- 
der non-resonant illumination. For each current-time trace, the time 
(tup), during which the current exceeds the threshold, is integrated 



and gives the values (Ltup) plotted in Fig. |2b,c. As shown in Fig. |2b, 
the resonance shows up at the photon energy detuning of 1 jieV at 
zero magnetic field, and starts to split into two diagonal arms with 
increasing magnetic field. It is due to the Zeeman eff'ect of one in- 
dividual Er^^ ion, with Ag ^ 4.8. Similarly, the Zeeman splitting 
of the resonance around 1,538.0 nm is studied as shown in Fig. |2f:, 
and the rectangular regions denoted by the darkest blue colour are 
not scanned. There appear to be two resonances with similar reso- 
nant wavelengths and the same g-factor diff'erence (Ag ^ 3.3) but 
with different signal intensity. This could be due to two individual 
Er^^ ions with the same site symmetry but with different capacitive 
coupling. Furthermore, the Zeeman splitting of the resonance around 
1,538.0 nm shows polarization dependence. As shown in Fig. |2t, the 
diagonal arm is weaker than the anti-diagonal one. By modifying the 
polarization of the light entering the cryostat, the diagonal arm was 
tuned to be stronger than the anti-diagonal one. The site symme- 
try of individual Er centres can be determined with the polarization 
dependence and a rotating magnetic field measurement. Spin selec- 
tive excitation even for degenerate spin states can be achieved with 
the maximum contrast of the polarization dependence, which allows 
spin readout without a magnetic field. 

The hyperfine structure is of great interest as the nuclear spin has 
long coherence times for quantum information storage^^^^^^. In ad- 
dition, it is a strong evidence for distinguishing between diff'erent 
ions as well as other defects. Erbium has six stable isotopes, among 
which only ^^^Er has a nonzero nuclear spin of 7=7/2, leading to 
eight nuclear spin states. At high magnetic field, the hyperfine inter- 
action can be treated as a perturbation of the Zeeman effect^'', so each 
electron spin state will split into eight sublevels due to the hyperfine 
interaction (Fig. |3^). At low magnetic field, the hyperfine interaction 
is comparable to the Zeeman eff'ect, so the sublevels will mix^^. 

In order to investigate the hyperfine structure of i^-'Er^^ ions, we 
implanted ^^^Er and ^^^Er with zero nuclear spin as control ions. We 
first study the photoionization spectrum of an Er^+ control ion with 
zero nuclear spin. The integrated time (Ltup) is plotted as a function 
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FIG. 2. The Zeeman effect of individual Er^^ ions, a, Schematic di- 
agrams showing the Zeeman spUtting and optical transitions of Er^+ 
ions in silicon. The splitting of the "^113/2 and ^Ii5/2 states depends on 
the site symmetry of the Er centre, b, The Zeeman splitting scan of 
the Er resonance with the centred wavelength of 1,535.8 nm. Each 
pixel stands for a current- time trace recorded for 50 s. c, The Zee- 
man splitting scan of the Er resonance with the centred wavelength 
of 1,538.0 nm. 



of the photon energy detuning, as indicated by the blue dashed line in 
Fig. [3J3. The same spectral asymmetry as that in Fig. [T^ is observed. 
We attribute the asymmetry to the correlation between the Stark shift 
of the ^Ii5/2 - ^Ii3/2 resonance and a broadening of the Coulomb peak, 
both of which are sensitive to fluctuating electric fields in the chan- 
nel. The fluctuating field is attributed to the laser excitation of trap 
states in or near the channel. Since we directly observe the eff'ect on 
the Coulomb peak it is possible to remove part of this broadening of 
the peak. After applying this correction (details in Methods), a min- 
imum FWHM spectral width of 50 neV is observed as indicated by 
the red solid line in Fig. l3t. To significantly reduce the linewidth fur- 
ther it will be necessary to reduce the density of trap states. As well 
as the electric field induced shifts the line is expected to be broad- 
ened through magnetic interactions with ^^Si up to tens of neV and 
with paramagnetic centres in the device. It is expected from analogy 
with observations in Er^+:Y2Si05^^ that applying a large magnetic 
field will suppress this broadening mechanism. 

In the following, we show the hyperfine structure of one i67gj.3+ 
ion. As shown in Fig. [3j:, the photoionization spectrum taken at high 
magnetic field (B = 0.14 T), reveals eight resonant peaks with the 
photon energy diff'erence of about 0.2 jieV between each other. The 
high spectral resolution allows the nuclear spin readout with potential 
for single- shot readout and manipulating the nuclear spin states. As 
the addressability doesn't rely on a specific magnetic field, the pho- 
toionization spectra are measured at a series of magnetic fields, as 
shown in Fig. |3ji. The Zeeman shift is subtracted to show evolution 
of the hyperfine interaction. At high magnetic field, eight significant 
peaks are observed all through (0.08 T < 5 < 0.14 T), as the hyper- 
fine interaction can be treated as a perturbation of the Zeeman eff'ect. 



At low magnetic field, multiple resonances show up (-0.04 T < 5 < 
0.06 T), revealing the mixing of the hyperfine sublevels, since the 
hyperfine interaction is comparable to the Zeeman eff'ect. 

The eight significant peaks, representing the eight diff'erent nu- 
clear spin states of ^^^Er, demonstrate that the resonances are due to 
the i6Vgj.3+ JQj^ rather than other ions or defects. These eight hyper- 
fine peaks (Fig. [3j:) correspond to the allowed transitions (zlM/=0) 
preserving the nuclear spin states, but it's still a question whether 
they are due to the AMs=0 or AMs=±\ transitions. As shown in 
Fig. [3^, we attribute them to the AMs^O transitions for two rea- 
sons. First, the energy diff'erence between the two most distant hy- 
perfine peaks is only about 1.7 jieV (Fig. [3]:), which is much smaller 
than the typical splitting energy of the zIM^^+l transitions of Er^^ 
ions. The electron paramagnetic resonance measurements of ^^ ^Er^^ 
ions in crystals show a splitting energy (2AEl) of about 30 [leW^I^, 
which corresponds to the A Ms =+1 transitions. Second, a ninth peak 
shows up beyond the region between the two most distant peaks (at 
the photon energy detuning of -1.1 jieV in Fig. |3j:). The ninth peak 
is much weaker than the eight peaks but still recognizable, which we 
attribute to a not-allowed transition. The energy of the not-allowed 
transitions (zlM/=+l) of Er^+ ions can exceed the region between the 
two most distant peaks of the allowed transitions, only in the case 
of the A Ms =0 transitions. Consequently, the eight significant peaks 
are attributed to the AMs =0 transitions, and the splitting energy is 
expressed as \AEh - AEi\-\.l jieV. 

Hybrid optical/electrical access to single spins of individual ions 
in a nano-transistor has been demonstrated, which is applicable for 
other defects in solids. Specifically, with an Er-implanted SET the 
photoionization spectroscopy allows real-time observation of single 
optical excitation events avoiding the bottleneck of photon collec- 
tion. Furthermore, high-resolution optical frequency addressing cir- 
cumvents the limitations due to thermal broadening in earlier electri- 
cal detection of impurity spins ^ Our findings open the way to opti- 
cally address and manipulate the electron and nuclear spin states of 
an individual defect in a solid beyond the nitrogen- vacancy centre in 
diamond. In addition, this hybrid optical/electrical technique boosts 
the micro structural study of ions in a semiconductor to a single- site 
level, including microscopic aspects, electrical and optical activity, 
etc. 

An approach that combines dopant ions (e.g. Er, P) with quan- 
tum optical control and semiconductor fabrication technologies rep- 
resents an attractive platform to realize a scalable quantum compu- 
tation and communication architecture. Such a system could consist 
of individual ions inside a ring cavity coupled with each other via 
photons, and nearby charge- sensing devices used to read out the spin 
states of individual ions and to control the coupling between ions by 
Stark tuning. The ring cavities can be connected by optical waveg- 
uides, which enable quantum information transfer between individ- 
ual ions in diff'erent ring cavities. Here we demonstrated the first step 
towards such a system, i.e. optical addressing of individual ions, and 
further improvement can be made by reducing the observed linewidth 
as discussed previously. However, there are essential questions to be 
addressed in the future, such as electron and nuclear spin coherence 
times of Er (P) ions, the influence of photoionization on nuclear spin 
coherence, and spin-photon entanglement. 

METHOD SUMMARY 

The devices are fabricated with the same technique as the previous 
stud3p5 After complete device fabrication, an Er:0 co-implantation 
(dose ratio 1:6 ) is performed with the implantation energy of 400 
keV and 55 keV, respectively. There should be approximately 30-40 
Er ions in the sensitive region of one Coulomb peak. Under the er- 
bium implantation conditions we used, the beam is estimated to have 
been composed of 70-80% ^^^Er and 20-30% ^^^Er. The devices are 
then annealed at 700° C in N2 for 10 minutes to remove the implanta- 
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FIG. 3. The hyperfine structure of an individual Er^^ ion. a, Schematic diagrams showing the hyperfine spUtting and the AMs=0 transitions 
(orange dashed Une) of ^^'^Er^^ ions at high magnetic field, and the optical transition of Er^^ ions with zero nuclear spin at zero magnetic field 
for comparison (black solid line), b, The photoionization spectrum of an Er^+ control ion with zero nuclear spin. The red solid (blue dashed) 
line represents the data with (without) removing the broadening, c, The photoionization spectrum of a single i6V£j.3+ ^^^ j^iq eight significant 
peaks correspond to optical transitions of the eight nuclear spin states of ^^'^Er (7=7/2). d, The contour plot of the photoionization spectra of 
the i6^Er^+ ion showing evolution of the hyperfine interaction. 



tion damage and to initiate the formation of Er centres. All the mea- 
surements are carried out in a liquid helium cryostat at 4.2 K. The 
laser beam, with 4-5 mW optical power, goes through a single-mode 
fiber and irradiates the sample with a diameter of about 1 mm. In the 
initial phase of the experiments (Figs. [T]and|2]), a commercial tunable 
laser with an external cavity is used. To keep a high precision, we 
set one centred wavelength with the motor- actuator, and sweep the 
wavelength around the centred wavelength with the piezo-actuator. 
In the high-resolution experiments (Fig. l3]), the wavelength of an- 
other laser is stabilized to about 0.01 pm, and a wavelength meter is 
used to compensate the thermal drift. 
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METHODS 

Details of the devices. The devices used in this study are n-p- 
n field-effect transistors with a polycrystaUine siUcon gate wrapped 
around the p-type siUcon channel separated by the gate dielectric. 
The p-type channel has a boron doping of 10^^ cm~^. After com- 
plete device fabrication, an Er:0 co-implantation is performed with 
the implantation energy of 400 keV and 55 keV and the ion fluence 
of 4 X 10^^ cm~^ and 3 x 10^^ cm~^, respectively. This leads to an 
Er:0 dose ratio of about 1:6 in the channel region. Under the er- 
bium implantation conditions we used, the beam is estimated to have 
been composed of 70-80% ^^^Er and 20-30% ^^^Er. The presence 
of both oxygen impurities ^^ and boron impurities^ ^ is known to en- 
hance the luminescence of the Er^^ ions in silicon. The 700°C post- 
implantation anneal is within the thermal processing window for Er 
centre activation in silicon^ ^ 

In the experiments, the device is biased below the threshold volt- 
age, and only the corner regions of the silicon channel go into 
inversioiW A peak of the I-Vg curve is due to the Coulomb block- 
ade in one of the two corner regions, where the current flows. The 
sensitive region is defined as the region, in which one elementary 
charge change can be detected, with taking the current noise and the 
transconductance of the Coulomb peak into account. The sensitive 
region of one Coulomb peak is estimated to be the corresponding 
channel-corner region with a dimension of 100 x 50 x 20 nm (length 
X width X height) for the device shown in Fig. la. Simulations of 
the ion implantation based on SRIM^^ show that there should be ap- 
proximately 30-40 Er ions in the sensitive region of one Coulomb 
peak. 

Experimental details and data analysis. All the measurements 
are carried out in a liquid helium cryostat at 4.2 K. The laser beam, 
with 4-5 mW optical power, goes through a single-mode fiber and 
irradiates the sample with a diameter of about 1 mm. In the ini- 
tial phase of the experiments (Figs. [T]and[2]), a commercial tunable 
laser with an external cavity is used. To keep a high precision, we 



set one centred wavelength with the motor- actuator, and sweep the 
wavelength around the centred wavelength with the piezo- actuator. 
The current- time traces in Fig. [TJ: and Fig. [T}i taken at two diff'erent 
photon energies are consistent with the photoionization spectrum as 
indicated by the green and red diamonds in Fig. [T^, respectively. For 
instance, the current mainly stays at the background level (0.4 nA) 
at the photon energy detuning of -5 jieV, while the current jumps be- 
tween two levels (1.8 nA and 0.4 nA) at the photon energy detuning 
of 4 |ieV. It is worth to note that the two-level trace (Fig. [TJi) sug- 
gests only one single Er^+ ion is ionized. Multiple ions with diff'erent 
capacitive coupling will lead to a current-time trace with more than 
two levels, while two ions with the same capacitive coupling will 
lead to a current-time trace with three levels once they are simultane- 
ously ionized. We attribute the charge displacement to the ionization 
of an Er^+ ion rather than the charge fluctuations of the trap states, 
based on the observation that all the Er^^ ions that we observed con- 
tribute to a shift of the Coulomb peak towards lower gate voltage. In 
the high-resolution experiments (Fig. l3]), the wavelength of another 
laser is stabilized to about 0.01 pm, and a wavelength meter is used 
to compensate the thermal drift. The asymmetry as well as part of 
the broadening of the resonant peak is removed by adding an photon 
energy off'set to the data, and then the time, during which the current 
exceeds the threshold, is integrated and gives the values plotted as 
the red solid line in Fig. [3]d,c. In comparison, the red solid line with 
removing the broadening shows smaller widths and less noise than 
the blue dashed line without removing the broadening, nevertheless, 
the resonances in the latter are still clearly visible, as shown in Fig. 
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